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Abstract 

We report results from an experiment measuring the semi-inclusive reaction D(e, e'ps) where the 
proton Ps is moving at a large angle relative to the momentum transfer. If we assume that the proton 
was a spectator to the reaction taking place on the neutron in deuterium, the initial state of that 
neutron can be inferred. This method, known as spectator tagging, can be used to study electron 
scattering from high-momentum (off-shell) neutrons in deuterium. The data were taken with a 
5.765 GeV electron beam on a deuterium target in Jefferson Laboratory's Hall B, using the CLAS 
detector. A reduced cross section was extracted for different values of final-state missing mass W*, 
backward proton momentum ps and momentum transfer Q^. The data are compared to a simple 
PWIA spectator model. A strong enhancement in the data observed at transverse kinematics is not 
reproduced by the PWIA model. This enhancement can likely be associated with the contribution 
of final state interactions (FSI) that were not incorporated into the model. A "bound neutron 
structure function" was extracted as a function of W* and the scaling variable x* at extreme 
backward kinematics, where effects of FSI appear to be smaller. For ps > 400 MeV/c, where the 
neutron is far off-shell, the model overestimates the value of in the region of x* between 0.25 
and 0.6. A modification of the bound neutron structure function is one of possible effects that can 
cause the observed deviation. 

PACS numbers: 24.85. +p, 25.30.-c, 21.45. +v 

Keywords: deuterium, off-shell, neutron, structure functions 
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I. INTRODUCTION 



Decades before the nucleon substructure was discovered, numerous models were devel- 
oped that successfully describe most nuclear phenomena only in terms of nucleons, their 
excited states and strong force mediators - mesons. Nucleons and mesons are often called 
the "conventional" degrees of freedom of nuclear physics. The fundamental theory of strong 
interactions, quantum chromodynamics (QCD), describes physical processes in terms of 
quarks and gluons. QCD is very successful in describing the interaction of quarks at short 
distances, where perturbative methods, similar to those of quantum electrodynamics (QED) 
in atomic physics, are applicable. However, the same perturbative methods cannot be ap- 
plied to solve QCD at the length scales of a nucleus. The present difficulty to make rigorous 
predictions based on QCD at low momenta (corresponding to large distance scales) leaves us 
no choice but to continue to employ nuclear theories based on "effective" degrees of freedom 
- nucleons and mesons. In an attempt to resolve this discontinuity of theories, the focus of 
modern nuclear physics has turned to the intermediate region where QCD is not yet solvable, 
but the quark-gluon substructure of the nucleons must be taken into account in the nuclear 
models. 

One example of the interface between a hadronic and a quark-based description is the 
(possible) modification of the (quark-) structure of a nucleon that is part of a tightly bound 
pair. Due to the Heisenberg uncertainty principle, large momenta of the nucleons inside 
the nucleus can be associated with small internucleon spatial separations. The kinematical 
conditions are particularly clean in the case of the deuteron, where the relative motion of 
the two nucleons is completely described by the wave function in momentum space, ip{p). In 
all models of the deuterium nucleus, the nucleons have mostly low momenta and therefore 
are relatively far apart. However, even in the wave functions obtained from non-relativistic 
models of the nucleon-nucleon potential, there is a probability for the nucleons to have 
momenta high enough so that the proton and neutron can come very close together or even 
overlap. In such high density configurations the quark distribution within a nucleon can 
become modified either through off-shell effects jl| or through direct modification of the 
shape and size of the nucleon It is also possible that under these conditions the 

nucleons start to exchange quarks with each other or even merge into a single "six-quark 
bag" 0,0]. The quark-gluon degrees of freedom thus might play a direct role in modifying 
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nucleoli structure in high-density nuclear configurations. The analysis presented here is 
aimed at advancing the understanding of high density, high momentum nuclear matter. 

To study these high density configurations, we have used electron scattering from a high- 
momentum nucleon within a nucleus. In the case of a deuteron target this can be easily 
verified by taking advantage of the inherently simple structure of the two-nucleon system. 
If all the momentum and energy is transferred to the neutron, the proton is a spectator 
to the reaction and recoils with its initial momentum. Assuming that the detected proton 
was indeed a spectator to the reaction, the initial momentum of the struck neutron can be 
obtained using momentum conservation. Thus the neutron is "tagged" by the backward going 
spectator proton (for a extensive discussion of the spectator picture see, e.g.. the papers by 
Simula p| and Meltnitchouk et al. [ll). Measurement of a high-momentum proton emitted 
backwards relative to the momentum transfer direction allows us to infer that the electron 
interacted with a high-momentum neutron in deuterium. 



II. THEORETICAL MODELS 

A. Nucleons in the Nuclear Medium 



Energy conservation applied to the deuterium nucleus requires that the total energy of 
the proton and neutron bound within a deuteron equals the mass of the deuterium nucleus: 

Ep + E„ = Ma. (1) 

At the same time, the mass of the deuteron is less than the mass of a free proton plus the 
mass of a free neutron, Ma = Mp + M„ — 2.2246 MeV. Therefore, both the bound neutron 
and proton can not be on the mass shell at the same time. In the "instant form" dynamics, 
one of the nucleons is assumed to be on-shell, while the other one is off-shell and its off-shell 



energy is E*^ = - ^ + pi 

The final state motion of the on-shell ("spectator") nucleon can be described by its mo- 
mentum Ps or the light cone fraction a^: 

where = {Es,Pt,Psu) is the spectator proton momentum 4- vector. The component Ps,, of 
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the proton momentum is in the direction of the momentum transfer g, and px is transverse 
to q. 

Using a non-relativistic wave function ipNuips), the "target density" of neutrons which 
are correlated with spectator protons of momentum Ps can be expressed as: 

P{Ps) = J-\i^NR{Ps)\\ (3) 

where J = 1 + ^ = ^^7""-'^°, is a flux factor that accounts for the motion of the struck 
nucleon. 

The probability P{Ps) is related to the spectral function: 

dcy. 

Sias,PT) — -(fpT = PiPs)d^Ps, (4) 
as 

which yields S = Eg ■ P{ps)- 

In the light-cone dynamics framework, a non-relativistic deuterium wave function can be 
rescaled to account for relativistic effects at high momenta j^: 

^^^(«s,Pr)— ciV = \^NRm'')?d?k (5) 
as = l- , ^" (6) 

M2 + fc2 



VT = kT \h\ = ^^^-M\ (7) 

where is the light-cone fraction of the nucleus carried by the spectator nucleon and 



= {ko, kx, k\i) is its internal momentum, with k^ = y + k"^. The relativistic effect, in 
this picture, manifests itself in that the measured momentum of the nucleon Ps^^ is rescaled 
in the lab frame from the internal momentum k\\. The resulting deuterium momentum 
distribution is given by the spectral function: 



S'^^{as,pT) = \ ^ l^PNRim'. (8) 
2- as 

The spectral function is normalized to satisfy the relation: 

S'^^ias,PT)—d'pT = l. (9) 

as 

In the PWIA spectator approximation, the recoiling proton is on-shell at the moment of in- 
teraction and receives no energy or momentum transfer, so that its initial and flnal momenta 



in the lab are the same. The differential cross-section on a moving nucleon (with kinematics 
defined by the spectator variables as,PT) can then be calculated as: 



do- _ 47ra^ 



dx'dQ'^ 



ILZ L n _ 0,*^ J- M*^x*^y*^ 1-R 



(10) 



_2{1+R) + + Q2 

xF2{x*,as,PT,Q^) ■ S{a„pT)^d'^PT 

where S{as,PT)^d?PT is the probability to find a spectator with the given kinematics. In 
this expression, F2{x* ^agiPTiQ"^) is the "off-shell" structure function of the struck neutron 
and i? = — is the ratio between the longitudinal and transverse cross sections. The asterisk 
is used for variables that have been defined in a manifestly covariant way. For instance, the 
Bjorken scaling variable x = and the variable y = ^ that are valid for the scattering 
from a free nucleon at rest are replaced with their counterparts for the scattering on a moving 
neutron inside the deuteron: 

* _ Q ^ Q (W) 

^ ~ 2p%qf' ^ 2Mz/(2 - a,) ~ 2 - a, ^ ^ 

y = ~ y, 

where = (z/, q) is the momentum transfer 4-vector, = {E, 0, 0, E) is the momentum 
4- vector of the incident electron, p% = (M^ — Eg, —ps) is the momentum 4-vector of the 
off-shell neutron and Ma is the mass of the deuterium nucleus. In this approximation the 
struck nucleon is assumed to be on the energy shell, but off its mass shell. The mass of the 
free nucleon M is therefore replaced with the off-shell mass of the bound nucleon: 

M*^ = {M,-E,f-pl (12) 

The invariant mass of the final hadronic state in D(e, e'ps)X scattering can be expressed as: 



W*^ = (K + g^)' = M*2 + 2{Md -Es)v + 2p,|| |gl 



where it was assumed that ~ 2M. In the (Bjorken) limit of 1 the fraction in the 

brackets of the last term in equation (fT3|) takes the familiar form of the light-cone fraction 



of the nucleus carried by the spectator proton = "^J"" , yielding: 

_ _ q2 ^ 2Mz/ (2 - as) . (14) 

If one assumes that F2 is equal to its on-shell form, F2{x*,as,PT,Q'^) = F2^'^^ {x* , Q"^) , 
and integrates over the spectator kinematics, one obtains the usual convolution result for 
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the inclusive nuclear structure function F2a- In this picture the nucleus is built from free 
nucleons, i.e. the struck nucleon has the same quark distribution as a free nucleon. Any 
observed modification of the cross section from that of a collection of free nucleons is just 
due to the kinematic rescaling (Eqs. EJ because of the motion of the nucleons inside the 
nucleus. However, the difference in the x dependence of the inclusive deep inelastic cross 
section for free and bound nucleons observed by the European Muon Collaboration (known 

Q 

as the EMC-effect [2|), cannot be interpreted solely in terms of such a kinematic shift. A 
large number of models have been proposed to explain the EMC-effect. A good review of 
this subject is given by Sargsian et al. in Ref. j^. 

The most conservative approach assumes that any modification of the bound nucleon 
structure function is solely due to the fact that the struck nucleon is off its mass shell (E < 
M); for example see Ref. Other models invoke a change of the nucleon size and therefore 
a rescaling of the structure function with momentum transfer Q^, as in Ref. j^. Frankfurt 
and Strikman ^ link the modifications to the structure function with a suppression of small 
(point-like) valence configurations of a strongly bound nucleon. The most unconventional 
attempt to explain the EMC-effect is that of Carlson and Lassila |J,l5fl where nucleons inside 
of a nucleus in its high-density configuration are thought to merge and form multiquark 
states. For the case of deuterium, as much as 5% of the wave function would be in a 6-quark 
state in this model. The cross section for backward proton production is then expressed 
as a convolution of the distribution function for the valence quarks in a 6-quark cluster 
V^^^ and the fragmentation function for the 5-quark residuum into a backward proton, 
Dp/5q{z) oc (1 — z)'^, with z = a/{2 — x). 

Although all of these models can describe at least some aspects of the EMC-effect, they 
predict considerably different changes of the internal structure of deeply bound nucleons. 
These changes are masked in inclusive measurements, where one averages over all bound 
nucleons, most of which are below the Fermi surface. By selecting tightly bound nucleon 
pairs (with a fast backward going spectator as "tag"), our experiment can study these possible 
modifications more directly. 
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B. Final State Interactions 



The PWIA picture described above has to be modified to include the effect of final state 
interactions (FSI) and two-body currents (meson exchange currents). According to existing 
models (see below), there are kinematic regions where FSI are thought to be small, and 
other regions where FSI are enhanced. Reliable models of FSI exist for nucleon-nucleon 
rescattering Q|. In the resonant and deep inelastic region, the estimation of FSI is a lot 
more challenging. FSI can be modeled by replacing the spectral function in Eq. Uni with a 
distorted one: S^'^^as^pr) ■ 

Melnitchouk, Sargsian and Strikman [l| use the eD epn reaction as a first estimate 
of FSI in electron scattering from the deuteron. This calculation shows that for > 2 — x 
and pt close to zero FSI are small. In this model S^^^ is evaluated using a distorted wave 
impulse approximation (DWIA). According to this paper, FSI effects should not strongly 
depend on x, thus the ratios of the cross section for different ranges in x should be a good 
tool to look for the EMC-effect in the semi-inclusive eD ^ epX process. In the limit of 
large x, FSI become much more important for heavier nuclei, where rescattering hadrons 
produced in the elementary deep inelastic scattering (DIS) off the short-range correlation 
are dynamically enhanced. Therefore, deuterium targets, in the authors' opinion, provide 
the best way of studying the origin of the EMC effect. 

A more recent publication by Cioffi et al. Q| discusses backward proton production and 
FSI associated with DIS by evaluating S^^^ within a hadronization framework. The reinter- 
action of the backward-going spectator protons with the debris formed in a hadronization 
process is modeled using an effective cross section: 

a^^ = a^N ^a-^^UM + nal (15) 

where a^^ and a'^^ are the total nucleon-nucleon and meson-nucleon cross sections, respec- 
tively, and um and are the effective numbers of created mesons and radiated gluons. The 
cross section asymptotically tends to exhibit a simple logarithmic behavior. The magnitude 
of the effective reinteraction cross section differs significantly for different models, especially 
at angles of proton emission ~ 90°. This kinematic region is proposed by the authors 
as the best place to test various models of hadronization. In contrast with the calculation 



discussed in the beginning of the section, the model of 



id] predicts significant FSI for proton 



momenta \ps\ > 250 MeV/c even at extreme backward angles. 
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III. EXISTING DATA OVERVIEW 



Few data exist on the semi-inclusive scattering of a lepton from deuterium with a recoil- 
ing nucleon in the backward direction with respect to the momentum transfer. The data 
published so far were taken using either neutrino or antineutrino beams and had very low 
statistics that do not allow detailed investigation of the cross sections of interest. These ex- 
periments (see Berge and Efremenko focused on measuring the momentum, energy, 
and angular distributions of protons in the backward hemisphere relative to the beam line. 
Despite the low statistics, a notable difference in the distributions for backward and forward 
protons was observed. The data were shown to agree well with a pair-correlation model in 
which the detected backward proton is assumed to be a spectator to the reaction. 

The cross section ratio a^^ ja^ measured by the European Muon Collaboration [3| (where 
a^^ and are cross sections per nucleon for iron and deuterium respectively) showed 
deviations from unity (now known as the EMC-effect) that could not be explained only 
in terms of nucleon Fermi motion. That was the first evidence that the nuclear medium 
influences DIS processes. It provided an indication that nuclear matter is getting modified 
as its density increases. The effect was later confirmed by data from SLAC ll^ and 
CERN Q. 

An independent measurement of the modification of the quark structure of nuclei was later 
done at Fermilab 3| using continuum dimuon production in high-energy hadron collisions. 



known as the Drell-Yan process [131 . The measurement has shown no nuclear dependence 
in the production of the dimuon pairs in the region 0.1 < x < 0.3, and therefore, no 
modification of the antiquark sea in this range. A number of models developed to explain 
the EMC-effect in terms of strong enhancement of the pion cloud were ruled out by this 
experiment. 

A .ecen. p„,a..a.i„„ t.«„.fe. by D.eteric. and St.auc. Q H H Q 

in the ^He(e, e'p)^H reaction suggested medium modification of the electromagnetic form 
factors of the nucleon. The observed 10% deviation from unity could only be explained by 
supplementing the conventional nuclear description with effects due to medium modification 
of the nucleon as calculated by the QMC model [2^ [2^. 

A model in which the neutron and proton form a single 6-quark cluster was recently tested 
0| against old backward proton production data from neutrino scattering on deuterium col- 
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FIG. 1: (Color online) CLAS event with forward electron detected in coincidence with a backward 
proton. 

lected at Fermilab [^l- These data had sufficient acceptance for backward protons but 
were not previously analyzed for this signal. The proton spectrum from neutrino and an- 
tineutrino scattering from deuterium, taken at CERN 2^, was also discussed. The authors 
compared the momentum distribution of backward protons with the prediction of a 6-quark 
cluster model. Predictions of the model were shown to be in good agreement with the data, 
however, the statistics of the data were not sufficient to study the dependence on any other 
kinematic variables. 

In summary, existing data on inelastic scattering off nuclei average over at least some 
of the relevant kinematic variables (x, Q^, and the momentum of the struck nucleon) and 
are often limited in statistics. Only a more detailed analysis of the dependence of the cross 
section on these variables can yield clear distinctions between different models and theoretical 
descriptions of nucleons bound in nuclei. The experiment on the reaction D(e, e'ps) described 
here is the first to collect sufficient statistics for this purpose. 



IV. EXPERIMENTAL SETUP 



The data were collected over a period of 46 calendar days in February and March of 
2002 at the Thomas Jefferson National Accelerator Facility (TJNAF). We used a 5.75 GeV 
electron beam with an average current of 6 — 9 nA. The experiment was staged in Hall B 
of the TJNAF, where the CEBAF Large Acceptance Spectrometer (CLAS) is installed. Six 
superconducting magnetic coils divide CLAS into six sectors symmetrically located around 
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the beamline. Each sector covers almost 60° in azimuthal angle and between 10° and 140° 
in polar angle, thus providing almost in acceptance for charged particles. CLAS sectors 
are equipped with identical sets of detector systems (Fig. H}: 1) three regions of drift 
chambers (DC) track charged particle's passage though the region of magnetic field; 2) a 
layer of scintillating paddles form the CLAS time-of-fiight system (TOF); 3) the Cherenkov 
counters (CC) are installed in the forward region (10° < 9iab < 50°) of the detector and 
efficiently discriminate electrons from pions up to the particle momenta p ^ 2.7 GeV/c; 4) 
several layers of lead and scintillating paddles form the electromagnetic calorimeter (EC) 
designed to separate electrons from minimum ionizing particles. CLAS is described in detail 
in Ref. 3- 

A conical cryogenic 5 cm target, installed in the center of CLAS, was filled with liquid 
deuterium at a temperature of 22 K and pressure of 1315 mbar with a density of 0.162 g/cm^. 
The average beam current of 8 nA produced a luminosity of 1.1 x lO^'' cm~^ ■ s"^. 

The CLAS trigger was formed by a coincidence between CC and EC. The signal level 
for the trigger coincidence was set to be at least 1 photoelectron in CC and 0.5 GeV in EC. 
The level 2 trigger required a DC track candidate in the sector of the calorimeter hit. With 
this trigger configuration, the data rate was about 3 kHz and the dead time was usually less 
than 13%. 

Out of 4.5 billion events collected over the experimental run, only 350 thousand contain 
an electron in coincidence with a backward proton. The typical event of that type detected 
in CLAS is shown in Fig. [H The collected data sample has wide coverage in kinematics 
of the electron and proton (Fig. The momentum transfer ranges between 1.2 and 
5.5 GeV^/c^, while the invariant mass covers the quasi-elastic, resonant and deep inelastic 
regions. Protons were detected at large angles relative to the momentum transfer vector g, 
up to angles of 9pg ^ 145° and with momenta above 0.28 GeV/c. 

V. DATA ANALYSIS 

In this section we discuss all the key analysis steps that led to the extraction of the final 
results. 
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FIG. 2: Kinematic coverage for electrons {W vs. Q^) (a) and for recoiling protons (momentum ppr 
vs. polar angle 9pg) (b), within fiducial cuts. 

A. Event Selection 

The focus of this analysis is the ed e'psX reaction, therefore events containing coinci- 
dences between the scattered electron and recoiling proton have to be selected first. 

The scattered relativistic electron is expected to be the first particle that arrives at the 
detectors after interacting with the target nucleus. The particle was identified as an electron 
if it was the first in the event and its charge was measured by the DC to be negative. Electron 
identification (ID) cuts on the response of two of the remaining detector systems, CC and 
EC, reduce the background of vr" in the electron spectrum. The CC are very efficient 
in pion rejection up to P ~ 2.7 GeV/c, where pions start to emit Cherenkov light. For 
lower momenta of the particle P < 3.0 GeV/c a software cut of 2.5 photoelectrons was 
required to identify an electron. For the part of the data with particle momentum P > 3.0 
GeV/c, a software cut of 1 photoelectron was used (and the fiducial region increased - see 
below) to increase acceptance. The electron produces an electromagnetic shower in the EC 
immediately after it enters, while pions make mostly a minimum ionizing signal with a small 
sampling fraction (E/P). The minimum ionizing particles can be easily rejected by requiring 
that the visible energy deposited in the first 15 layers of the EC is ECinner > 0.08 ■ P and 
the total visible energy in the EC is ECtotai > 0.22 ■ P. 

In order to reduce the systematic uncertainty in the quality of electron identification, 
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detector fiducial cuts are applied. The fiducial region of CC is known to be within the limits 
of the EC fiducial region; therefore only a CC cut needs to be applied. We defined the 
fiducial region such that the CC was at least 90% efficient. 

In addition to the particle charge information, the DCs also measure the length from 
the target to the TOF system and the curvature of the track. From the curvature of the 
track the particle momentum can be reconstructed. The proton is identified using TOF time 
measurement (^tof) and DC momentum {poc) and track length (r) information. Assuming 
a positively charged particle is a proton, its velocity is given by 



where Mp is proton mass. Then the time the proton travels from the target to the TOF is 
toe — f/vDC- The particle is identified as a proton if the time difference At = toe ~ troF, 
corrected for the event start time, is within a time window —2 ns to 7 ns. 

A vertex cut is applied to ensure that the interaction took place within the volume of the 
target. The electron was required to have a vertex —2 cm < Z^i < 1.5 cm while the proton 
vertex cut was set to —2.5 cm < Zp^. < 2 cm (the target extends from -2.5 cm to 2.5 cm). 
Additionally the vertex difference between Zei and Zpr was required to be less than 1.4 cm 
to reduce the background from accidental coincidences. 

B. Kinematic Corrections 

The geometrical and structural complexity of CLAS is responsible for minor discrepancies 
in the measurement of the momentum and direction of a particle. These discrepancies are 
thought to be primarily due to the uncertainty in the magnetic field map and DC position. 
The effect of a displacement of the drift chambers and possible discrepancies in the measured 
magnetic field on the measured scattering angle 6rec and momentum p can be parameterized. 

The correction function contains 8 parameters describing the drift chamber displacements 
and rotations and 8 parameters describing the possible uncertainties in the magnitude of 
the magnetic field on the path of the particle. These parameters can be determined using 
multi-particle exclusive reactions which are fully contained within the CLAS acceptance. In 
an exclusive reaction all of the products of the reactions are detected and no mass is missing. 
Therefore, the kinematics of the reaction are fully defined and the goodness of fit can be 



Pdc 



(16) 




15 



evaluated using momentum and energy conservation. More details on this method can be 
found in Ref. (221 . 

For low-energy protons (P < 0.75 GeV/c) energy loss in the target and detector is signif- 
icant and needs to be corrected for. This energy loss was studied with the CLAS GEANT 
simulation and an appropriate correction was applied to the data. 



C. Backgrounds 

Even after the ID cuts described above, pions remain a non-negligible background in the 
electron spectrum. Their contribution needs to be estimated and appropriate corrections 
applied to the data. This was done using a sample of pions within EC cuts of Einner < 0.05 
GeV and Etotai < 0.1 GeV. The spectrum of photoelectrons in the Cherenkov Counters of this 
pion sample was scaled such that the sum of the normalized spectrum and that of a "perfect" 
electron sample (from a simulation normalized to data within a tight EC cut) agreed with 
the measured Cherenkov spectrum for electron candidates within our regular EC cuts. This 
normalized pion spectrum was then integrated above the software ID cuts of 2.5 and 1.0 
photoelectrons (depending on the data momentum range) and used to estimate the fraction 
of pions remaining in the electron sample after the Cherenkov ID cut. This fraction was 
fit to an exponential in pion energy and the resulting estimate of the pion contamination 
(ranging to no more than 6%) was used to correct the extracted data. 

A similar technique was used to measure the rate of positrons relative to that of electrons, 
by taking positive charge tracks and fitting their energy spectrum in the EC with a com- 
bination of "pure" pions (based on Cherenkov response) and "golden electrons" (very high 
Cherenkov cut). This positron to electron ratio can be used to estimate the fraction of the 
detected electrons which were not scattered from the beam but came from pair production 
7 e~^e~ or the Dalitz decay 7r° — 'ye~^e~. Once again, an exponential fit to the ratio 
was used to estimate this contamination for all kinematic bins and correct our final data 
accordingly. 

Despite the vertex cuts there is still a chance of having an accidental coincidence between 
an electron and a proton in the data sample. The background of accidentals has to be 
estimated and subtracted. At the same time, the loss of "true" protons due to the time and 
vertex cuts has to be determined. A purely accidental proton was defined as a positively 
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charged particle with the time-of-flight measured by the TOF to be at least 12 ns longer 
than the expected time-of-flight of a proton with that momentum. The time window for the 
accidental proton was taken to be 9 ns, the same as the proton ID time window, so that the 
expected arrival time for the accidental proton would not be more than 21 ns different from 
the expected arrival time of the real proton. In the case where the time window of accidentals 
is less than 5 ns away from when the deuteron (from elastic scattering events) would have 
arrived at the TOF counter, the accidental proton is deflned to be within a 9 ns window 
starting at 5 ns after the expected arrival time of a deuterium ion. The average background 
of accidental coincidences per nanosecond of the proton time vertex was calculated from 
the rate in the "accidental time window" described above and compared with the unbiased 
data sample of coincidences with good proton PID. The level of understanding of accidentals 
was also tested using the simulation results. The sum of the measured accidentals and the 
simulation is in agreement with the data on good electron-proton coincidences as selected 
by PID cuts (Fig. IS]). A small discrepancy on the positive side of the AZ distribution is due 
to another type of unwanted coincidences where a particle originating from the flrst electron 
vertex reinteracts further along the target cell, liberating a (backward) proton which arrives 
on-time with respect to the TOF. Protons produced in such a way enhance the positive 
side of the vertex difference distribution. The selected sample of accidentals contains only 
off-time events, and therefore does not fully reproduce the shape of the vertex difference 
distribution. A properly scaled sample of these excess events was added to the sample of 
purely accidental coincidences deflned using off-time protons. 

D. Simulation 

To extract absolute results from our experimental data, the detector acceptance has to be 
evaluated and an appropriate correction applied to the data. An idealized model of all the 
detector systems of CLAS is implemented in the code known as "GSIM". The program is built 
on the foundation of the GEANT simulation software package, supported by CERN. GSIM 
allows simulation of the detector response to a propagating particle, simulating energy loss 
as well as emission of secondary particles during the passage of the particle through parts 
of the detector. After the response of the ideal detector is simulated, existing detector 
inefflciencies are introduced. This is done using a separate program called "GPP" (GSIM 
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FIG. 3: (Color online) Data for the difference between the electron and proton vertex (triangles) 
compared to a fit (solid histogram) composed of a simulation of true coincidences (not shown) and 
measured accidental coincidences (dash dotted histogram). The vertical dashed lines indicate the 
cut used to select data for analysis. 

post-processor). GPP uses precompiled information on dead regions of the DC and TOF to 
remove the signal for these parts of CLAS from the GSIM output. The final output is then 
analyzed exactly the same way as the real data. 

The events used as input for the CLAS GSIM simulation were generated follovi^ing the 
cross section Eq. Uni The Paris v^^ave function 0| v^as used to select the momentum of the 
"spectator" nucleon first. A comparison vv^ith the Argonne V18 vi^ave function j^l shovi^ed a 
negligible difference in the momentum distributions. The generated nucleon momentum can 
either be directly used follov^^ing the prescription for the non-relativistic spectral function 
(Eqs. IH|4|1 or as the "internal momentum" in the light cone description, Eqs. EHHl From the 
spectator nucleon kinematics, we then calculate the initial four-momentum of the struck 
nucleon and determine the scattered electron kinematics in the rest frame of that nucleon, 
then transform it back to the lab frame. That way, all of the "starred" variables in Eq. Uni 
are automatically evaluated with the proper relativistic rescaling. 

The electron scattering cross section used to generate the electron kinematics is based on 
the code RCSLACPOL that was developed at SLAC j^^l- It uses parametrizations of world 
data on unpolarized structure functions and elastic form factors. These parametrizations are 
described in l31 



and SLAC 



and are based on fits to unpolarized structure function data from NMC [s^l 
0, 0|- The nucleon form factors were taken from Ref. 01 • All form 
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factors and structure functions for bound nucleons are assumed to be equal to the free ones 
at the corresponding values of x (in the DIS region) or W (in the resonance region, with 
a smooth transition between both). The free neutron structure function was extracted 
from fits to the world data on the deuteron in a self-consistent manner by ensuring that our 
model, integrated over all spectator kinematics and summed over both proton and neutron 
contributions to electron scattering, agrees with those fits. 

Three different versions of the code were compiled to satisfy our needs for simulation of 
electron scattering on ^H: 1) elastic scattering on one nucleon in the deuteron (with the other 
being a spectator), including the elastic radiative tail; 2) inelastic scattering on one nucleon 
in the deuteron (with and without radiative corrections); and 3) elastic scattering off the 
deuteron nucleus as a whole. Radiative effects can be included in the simulation following 
the p.esc.ip.io„ MO and Tsai Q. In the fl.. Wo ca.es, these radiative co.eetions a.e 
applied to the electron scattering cross section for the struck nucleon in its rest frame, while 
the spectator simply determines the kinematic transformation into the lab system. The 
generator is capable of simulating both inclusive D(e,e') (by adding the first two processes 
for both protons and neutrons with the third one) and semi-inclusive D(e, e'p^) processes, 
which is controlled by a configuration file. While this generator may not be very realistic 
in its description of the underlying physical processes (since it does not contain FSI, non- 
nucleonic currents in deuterium, or modifications of the nucleon structure function for off- 
shell nucleons), it is sufficiently accurate (see below) to allow a largely unbiased extraction 
of the acceptance and efficiency of CLAS, by comparing accepted simulated events to the 
initial distribution of generated events. 

The quality of the simulation procedures can be evaluated by comparing the predicted 
number of counts for well-studied processes in data and simulation. To date, one of the 
best studied cross sections in nuclear physics is that of elastic electron scattering from a free 
proton. To select elastic events a cut on the invariant mass W was used: 0.9 < W < 1.1 GeV. 
The overall shape is reproduced well and the measured cross section lies well within 10% 
of the simulated one at low (where our statistical error allows a significant comparison). 
The distribution of the simulated inclusive cross section for quasi-elastic scattering on 
deuterium is also in good agreement with the experimental data. Here the events were also 
selected using the invariant mass cut 0.9 < W < 1.1 GeV. In the region of relatively good 
statistics at low the deviation from unity on the data to simulation ratio does not exceed 
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10%. Finally, the rate of inclusive D(e,e')X events for all final state invariant masses W 
agrees with the prediction of our model to within 5-10%. 

A sample of simulated events that exceeds the statistics of the experimental data by a 
factor of 10 was generated for the D(e, e'ps) reaction and was used in the analysis to correct 
the data for detector acceptance and bin averaging effects. The high event count of the 
Monte Carlo assures that the statistical error of the data points are not dominated by the 
statistical error of the simulation. 

E. Result Extraction 

The events from the data set were sorted in four-dimensional kinematic bins in W* (or 
X*), Q^, Ps and cos9pg (or and pt). We chose two bins in Q^, one with 1.2 (GeV/c)^ < 
< 2.1 {GeV/cf (average = 1.8 {GeV/cf) and one with 2.1 {GeV/cf < < 5.0 
(GeV/c)^ (average = 2.8 (GeV/c)^, and five bins in ps, with average values of ps = 
0.3,0.34,0.39,0.46 and 0.53 GeV/c. 

To extract the final results, the above bins were filled separately for the following cate- 
gories of events: 1) experimental data with all the standard electron and proton ID cuts; 
2) accidental electron-proton coincidences based on experimental data; 3) coincidences with 
protons from secondary scattering events; 4) simulated data for the elastic scattering on a 
bound neutron, including the radiative elastic tail; 5) simulated data for the inelastic scatter- 
ing on a bound neutron. Accidental coincidences and coincidences with secondary protons 
were then subtracted from the data on a bin-by-bin basis. The simulated elastic scattering 
data were also used to subtract the elastic radiative tail from the experimental data. For 
this purpose both data and simulation were first integrated in the range of the invariant 
mass of the unobserved final state W* from 0.5 to 1.1 GeV. The elastic radiative tail in the 
simulation was then scaled by the ratio of the data to the simulation and subtracted. 

As was previously discussed, in the spectator picture, the cross section for the off-shell 
nucleon can be factorized as a product of the bound nucleon structure function and the 
nuclear spectral function, multiplied by a kinematic factor (see Eg. fTO|) . Using the data of 
this experiment, it is possible to extract this product, and, in the region where FSI are small 
and the spectral function is well described by the model, even the off-shell structure function 
by itself. To do that, the experimental data (with accidentals, rescattered proton events, 
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and elastic radiative tail subtracted) were first divided by the simulated inelastic data. 
The simulated events were generated using the cross section Eq. [101 with full consideration 
of radiative effects. To extract the product of structure and spectral functions, the ratio 
of data to simulation was multiplied with the product F2n{x*,Q'^) x S{as,pT), calculated 
using the same model that was used in the generator. Similarly, to obtain the product of 
the structure function with the probability distribution for the proton momentum in 
deuterium, we multiplied the ratio of data to simulation with the factor F2„(x*, Q"^) x P{ps) 
from our generator model. In both cases, the dependence of the extracted data on the specific 
model for the simulation is minimized, since the "input" {F2n and S{as,PT) or -P(ps)) cancels 
to first order. Basically, this procedure corrects the data for the detector acceptance, bin 
migration and radiative effects, and produces a "normalized cross section" by dividing out 
the kinematic factor '^'^*q¥ as well as the factor in square brackets in Eq. [101 (which depends 
weakly on the ratio R = ai/aT). To extract the ("off-shell") structure function the ratio 
of data to simulation was multiplied with the free nucleon structure function F2„(x*,(5^). 
This assumes that the spectral function used in the simulation describes the momentum 
distribution of the spectator protons reasonably well. 

F. Systematic Uncertainties 

To simplify the statistical error calculation, all the corrections for the detector inefficien- 
cies and data sample contamination (except for accidentals and the radiative elastic tail) 
were applied to the simulated events. 

The efficiency of the CC electron ID cut is well reproduced in the simulation. A 1% 
systematic uncertainty enters here to account for the observed deviation of the cut efficiency 
from sector to sector. The EC ID cut efficiency is reproduced only partially. The efficiency of 
the cut in data was found to be 95%, however the same cut, applied to the simulation, is 98% 
efficient. The difference might be a result of data being contaminated with pions, despite 
the increased CC threshold. The simulated data were scaled down by a constant factor of 
0.97 to account for the difference in the effect of the cut. A 2% systematic uncertainty was 
assigned to this factor due to the uncertainty about the source of the deviation. A variable 
factor that ranges from 1.06 to less than 1.01 was used to introduce pion contamination 
into the simulation. The factor varies with the particle scattering angle and momentum. 
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A variable factor was also applied to the electron spectrum in the simulation to introduce 
electrons coming from electron-positron pair creation. The resulting systematic uncertainty 
was estimated by varying these factors by 50% of their deviation from unity. The resulting 
change in the distribution in each of the final histograms was used as an estimate of the 
systematic uncertainty of these corrections. 

Some additional corrections were applied to the proton spectrum. A constant factor of 
0.99 was introduced to reflect the difference in the effect of the proton timing ID cut on 
the real versus the simulated data. The systematic uncertainty of 0.5% on this number 
accounts for the momentum dependence of the effect. A factor dependent on the proton 
momentum was applied to the simulated data to account for the discrepancy between data 
and simulation in the effect of the cut that was set on the difference between the electron 
and proton vertices. The systematic uncertainty here is evaluated individually for each 
histogram, by varying the correction by 50%. 

A major contribution to our systematic error comes from remaining differences between 
the simulated and the "true" inefficiencies of CLAS. Even after removing bad channels and 
accounting for all known detector problems, we find that the ratio of simulated to measured 
rates for reconstructed protons varies from sector to sector. We use the RMS variation 
between sectors to estimate this systematic error as about 11% on average. We also include 
a 3% scale error on the target density, effective target length, and beam charge calibration. 

The data were corrected for the radiative elastic tail and accidental coincidences by direct 
subtraction of normalized (simulated or real) data (see previous subsection). The normal- 
ization factors were varied by 50% of their deviation from unity to estimate the systematic 
errors due to these corrections. The uncertainty on the inelastic radiative corrections was 
also calculated as 50% of the deviation from unity of the correction factor. We checked our 
radiative correction procedure against the existing code "EXCLURAD" |2S| for the case of 
quasi-elastic scattering (pn final state) and found good agreement within the stated uncer- 
tainties. 

A final systematic uncertainty comes from the model dependence of our simulated data. 
While the model input cancels in our extracted values for F2n{x*,Q'^) x S^as^pr) to first 
order, both migration between adjacent kinematic bins and distribution of events within a 
bin (where the CLAS acceptance might vary) are somewhat model-dependent. We estimated 
this effect by modifying the model input to agree with the cross section extracted from our 



22 



Source of Uncertainty 


Typical Range (in % of data yalue) 


EC ID Cut 


2 


Trieffer Efficiency 


2 


Secondary Electrons 


0.7 


Electron Vertex ID Cut 


0.6 


Proton Timing ID Cut 


0.5 


CC Efficiency 


1 


Pion Contamination 


0.5 ... 3 


/e~ Contamination 


... 0.75 


Pure Accidental Coincidences 


... <1.2> ... 4 


Coincidences with Knock-out Proton 


... <2.3> ... 6 


Vertex Difference Cut 


0.75 ... 1.5 


Quasi-elastic Radiatiye Corrections 


... <1.9> ... 11 


Inelastic Radiative Effects 


... <2.7> ... 12 


Luminosity 


3 


Tracking Inefficiency 


11 


Bin Migration & Model-Dependence of 
Acceptance 


... <5.2> ... 10 


Total 


15.5 ... <16.9> ... 34.1 



TABLE I: Systematic errors in percent of the data values. The typical range of the error as well 
as their RMS values (in brackets) are given. 

data. The deviation of the simulated events with this modified cross section from the data 
is a direct measure of the magnitude of this systematic error. We found its magnitude to be 
generally below 5%, going up to 10% for higher proton momenta. 

All systematic errors were added in quadrature and are shown as shaded bands in the 
Figures in the following section. The summary of systematic uncertainties is presented in 
Table m 
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FIG. 4: (Color online) Data (points) and results of the Monte Carlo (MC) simulation based on two 
different PWIA models (solid and dashed curves) for the total number of counts versus cos 9pg for 
proton momenta ps = 280-320 MeV/c (a) and ps = 360-420 MeV/c (b), integrated over electron 
kinematics. The total systematic error is indicated by the shaded band. 

VI. RESULTS 

In the following, we show several representative histograms (one-dimensional projections 
of the four-dimensional bins), comparing our data to our simple PWIA spectator model to 
elucidate some general trends. 

In Fig. 131 we show as a first step the accumulated number of protons (in coincidence with 
a scattered electron) for several bins in cos 6pq, where 9pq is the angle between the virtual 
exchanged photon and the proton. The data are not corrected for acceptance and efficiency 
and therefore fall off at large angles where CLAS has limited acceptance. The curves shown 
are from our simulation of these data, including the CLAS acceptance and without any 
normalization. Using the light cone prescription (Eq. ^ for the momentum distribution 
of the initial proton (solid curve), good agreement between the data and our Monte Carlo 
(MC) simulation is observed up to cos 9pq ^ —0.3. The result for the non-relativistic wave 
function (Eq. El dashed line) is similar in these kinematics. At more forward angles the data 
exceed the simulation by a large factor, especially at higher momenta (Fig. IIJd), indicating 
a breakdown of the pure PWIA spectator picture. We assume that this enhancement is due 
to FSI between the struck neutron and the spectator proton (see below). 
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FIG. 5: (Color online) Momentum distribution of the recoiling proton. Data (points) are compared 
with our MC simulation (solid curve) for the range of recoil angle —1.0 < cosOpq < —0.3 (a) and 
—0.3 < cos^pg < 0.3 (b). All events within a missing mass range 1.1 < W* < 2.0 GeV were 
summed together for this plot. 

The momentum distribution plotted separately for backward {9pg > 108°) and transverse 
(72° < 9pg < 108°) proton kinematics confirms this picture for the relative importance of 
non-PWIA processes (Fig. E]). The momentum distribution of the backvv^ard protons is 
reasonably weW described by the PWIA model, indicating that distortions due to FSI are 
rather small in this region. At the same time, the momentum distribution for the transverse 
protons is strongly enhanced at momenta above 300 MeV/c, as predicted by several models 



of FSI 



Id 



40l |. For momenta below about 300 MeV/c, the acceptance and efficiency of 
CLAS for protons falls off even faster than predicted by our Monte Carlo simulation. This 
explains the fall-off at low momenta in Fig. [HI 

In Fig. El we look at the angular distribution of the protons in more detail. The reduced 
cross section described in the previous section is plotted for three different proton momenta 
(increasing from left to right), as well as three different missing mass ranges of the unobserved 
final state (increasing from top to bottom) in the reaction D{e,e'ps)X . Several trends can 
be observed: 

• At proton momenta around 300 MeV/c, the extracted reduced cross section is con- 
sistent with our simple PWIA spectator model throughout the whole angular range 
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FIG. 6: (Color online) Results for the normalized cross section (equivalent to the product x 
P{ps) in the spectator picture) for the reaction T){e,e' ps)X . Each row is for a different mass W* 
of the unobserved final state X, namely W* = 0.94 GeV (quasi-elastic scattering) in the first row, 
W* = 1.5 GeV in the second and W* = 2 GeV in the third. The three columns are for three different 
proton momentum ranges, with average momenta of ps = 0.3, 0.39 and 0.56 GeV/c, respectively. 
All data (filled circles with statistical error bars) are for our lower bin (with average of 1.8 
(GeV/c)^). The two lines come from our simple PWIA spectator model using a light-cone wave 
function (solid line) or a non-relativistic WF (dashed line), while the shaded band at the bottom 
indicates the systematic error. 



and for all final state masses. This is consistent with expectations that destructive 
and constructive interference effects between FSI and PWIA cancel roughly in this 
momentum range [fll l4fl|. 

• For larger proton momenta, deviations from PWIA behavior show up as an increase of 
the normalized cross section at transverse kinematics. This increase appears approx- 
imately around cos 9pq = —0.3 and continues beyond cos 9pq = {9pq = 90°). Such 
an increase is not likely due to uncertainties in the deuteron wave function, which is 
isotropic in the non-relativistic case and is equal to the non-relativistic wave function 
for transverse proton momenta if one uses light-cone wave functions. However, such 
an effect is expected within models of FSI due to the initial motion of the nucleon on 
which the rescattering occurs (see Fig. 3 in Ref. and Ref. [40]). The strength of 
FSI in these models is the largest for the highest recoiling proton momenta, consistent 
with the trend of the data. 

• The non-PWIA effects seem to be more pronounced for the largest missing masses (see 
also below). This behavior is in qualitative agreement with the FSI model by Cioffi 
delgi Atti and collaborators (jiJ, |41|, where the strength of rescattering is related to 
the number of hadrons in the final state (Eq fTBll . 

This last point can be seen more clearly in Fig. which shows the ratio between the 
observed cross section and the prediction of our PWIA spectator model for proton momenta 
around 0.46 GeV/c, for four different ranges in final-state missing mass (slightly offset from 
each other for each point in cos6pg). The data for different missing mass values are statis- 
tically close to each other (and close to unity) in the backward region where rescattering 
effects can be assumed to be small. Conversely, in transverse kinematics the ratio substan- 
tially exceeds one and is largest for the highest W* bin. The enhancement in transverse 
kinematics is also large in the A— resonance region. This could be due to A-production in 
FSI between the struck neutron and the "spectator" proton. 

Concluding that the spectator PWIA model works reasonably well in the region of large 
backward angles {cos9pq < —0.3), we concentrate on this region to study the momentum (off- 
shell) dependence of the effective electron scattering cross section on the bound neutron. At 
first, we directly compare the extracted effective structure function of the off-shell neutron, 
for inelastic final states {W* > 1.1 GeV) to the on-shell structure function (see Fig. |HI). 
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FIG. 7: Ratio of data to model as a function of cos^pg for four values of missing mass W* at 
Ps = 460 MeV/c and = 1.8 GeVVc^. 

To obtain this structure function, the measured cross section was divided by the proton 
momentum distribution, Mott cross section and the kinematic factor as explained in the 
previous section. Even within the PWIA picture, the results could have a pg- dependent 
scale error because our simple model may not describe the nucleon momentum distribution in 
deuterium perfectly; however, the x*-dependence in each individual panel would be largely 
unaffected by such a scale error. Indeed, the data agree reasonably well with the simple 
parameterization of the free neutron structure function from our model at the two lower 
momenta (with average deviations of ±10%). At the higher two momenta, the data fall below 
the model in the range of x between 0.3 and 0.6 by as much as 20% - 30%. Such a reduction 
in the structure function is expected in several models of modification of bound nucleon 
structure Some residual FSI might also contribute to the observed x*-dependence, for 
instance by enhancing the region of small x* (corresponding to large W*). 

To reduce the model dependence of such comparisons as in Fig. [HI the authors of Ref. [l| 
suggested to take the ratio between the extracted "off-shell" structure function at some 
relatively large value of x* (where most models predict the biggest off-shell effects) to that 
at a smaller value of x* where the EMC-effect is known to be small. This ratio (normalized 
to the same ratio for the free neutron structure function, F2„) is plotted in Fig.lHlfor a range 
of transverse momenta 0.25 GeV/c < pr < 0.35 GeV/c. Nearly all dependence on our model 
cancels in this ratio; only the overall scale depends on the ratio of for free neutrons at 
two different values of x, which is not perfectly well known. The ratio plotted in Fig. [HI is 
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FIG. 8: (Color online) Results for the extracted "off-shell" structure function of the neutron 
in the PWIA spectator picture. The model (solid curve) is a simple parameterization of the free 
on-shell neutron structure function, modified to account for the kinematic shift due to the motion of 
the "off-shell" neutron. The sections of the plot correspond to different recoiling proton momenta: 
Ps = 300 MeV/c (a), ps = 340 MeV/c (b), p^ = 460 MeV/c (c) and Ps = 560 MeV/c (d). The 



quantityp 
Simula j^l- 



)lotted here is similar (but not identical) to the quantity F^'^ * ) defined in the paper by 



also independent of the deuteron momentum distribution P{ps)] however, according to some 
models FSI effects could be different for different x* . This seems to be born out by 
Fig. ini While all PWIA models of off-shell effects predict unity for the ratio at values of the 
light cone variable around 1, we find a strong suppression in the region up to ~ 1.1 
(corresponding to Opq around 90°) where FSI are most pronounced. This behavior could be 
explained within the FSI model of Ref. [l^ which predicts larger FSI effects for final states 
with a larger number of hadrons, leading to an increase of the denominator (cross section 
at small x* , which corresponds to large energy transfer to the unobserved final state). 
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FIG. 9: (Color online) Ratio of the extracted "off-shell" structure function F2n at x = 0.55, = 2.8 
(GeV/c)^ to that at x = 0.25, = 1.8 (GeV/c)^, divided by the ratio of the free structure functions 
at these kinematic points. The error bars are statistical only and the shaded band indicates the 
overall systematic error. This plot is for similar (but not identical) kinematics as Fig. 6 in the 



paper by Melnitchouk et al. 



Beyond ~ 1.1, the data still lie below unity (by about 17%) but appear fairly constant 
with tts- Although this suppression could be interpreted as an off-shell effect, the data 
appear inconsistent with some of the more dramatic predictions of a steep falloff for the 
ratio at high {e.g., Ref. j3|). The prediction for this ratio from the 6-quark cluster 
model Pi varies between 0.7 and 1 at = 1.4 and is therefore compatible with our result. 
Once realistic calculations including FSI effects become available for the kinematics of our 
data set, a more quantitative comparison with various models for the off-shell behavior of 
the structure function F2{x*, Q'^,Ps) will be feasible. Such calculations are underway |4^I3|- 



VII. SUMMARY 



Taking advantage of the large solid angle acceptance of the CEBAF Large Acceptance 
Spectrometer, a large amount of data (~ 350K events) was collected on the reaction 
D{e,e'ps)X in the exotic region of extreme backward proton kinematics. The data range 
from 1.2 to 5 (GeV/c)^ in momentum transfer and reach values of the missing mass of 
the unobserved final state W* of up to 2.7 GeV. Protons with momentum Ps as low as 
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280 MeV/c and up to 700 MeV/c were detected, at angles 6pq relative to the direction of the 
momentum transfer extending up to more than 140°. In terms of the light cone variables, 
the data span values of the light-cone fraction as up to about 1.7, with a minimum proton 
transverse momentum relative to q of 150 MeV/c and up to 600 MeV/c. 

Reduced cross sections were extracted as a function of W* (or Bjorken- variable x*) and 
«T, Vt (or cos6'pg, ps), for two large bins in Q^, allowing us to test theoretical calculations 
against the presented data. Comparison with a simple PWIA spectator model shows mod- 
erately good agreement in the kinematic region of lower momenta and cosOpq < —0.3. For 
increasing "spectator" momenta ps > 0.3EG V/c FSI and other non-PWIA effects become 
strong, especially in the region of proton scattering angles cos 9pq > —0.3. These effects 
seem to depend on the invariant mass W*; on the other hand, no strong dependence of these 
effects on momentum transfer is observed. This behavior is in qualitative agreement with 
models Q, |4l| that describe the strength of FSI in terms of the number of hadrons in the 
final state X. The angular {9pq) and momentum (ps) dependence of the observed strength 
in the cross section in the quasi-elastic region (where X is a neutron in its ground state) are 
also in good agreement with detailed calculations [l^ showing a transition from destructive 
interference below Ps = 300 MeV/c to a strong enhancement at Ps > 400 MeV/c around 
cosOpq = 0.2 (see Fig. El and also Ref. Q)- 

A depletion compared to the PWIA model is observed in the data at cos 9pq < —0.3 and 
for high Ps, where the struck neutron is far off its mass shell. This reduction might be due 
to nucleon structure modifications. It is especially apparent in the region of moderate x* 
which overlaps in part with the nucleon resonance region. However, it is also possible that 
our simple model predicts too much strength in the deuteron momentum distribution at 
these higher momenta. This would lead to an "apparent" depletion for all values of x* (or 
W*), which would be somewhat modified by a remaining FSI-induced enhancement at high 
W*. 

Ultimately, our data will serve to constrain detailed theoretical calculations, including off- 
shell and FSI effects. Once these effects are well-understood at high spectator momenta, one 
can safely extract the neutron structure function at lower momenta where those corrections 
are smaller and where their uncertainty will not affect the result. This method will be used 
in the upcoming "BoNuS" experiment at Jefferson Lab. A statistically improved data set 
with much larger kinematic coverage can be obtained once Jefferson Lab has been upgraded 
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to 12 GeV beam energy. 



Acknowledgments 

We would like to acknowledge the outstanding effort of the Accelerator, Target Group, 
and Physics Division staff at TJNAF that made this experiment possible. 

This work was supported by the U.S. Department of Energy, the Italian Istituto Nazionale 
di Fisica Nucleare, the U.S. National Science Foundation, the French Commissariat a 
I'Energie Atomique, the French Centre National de la Recherche Scientifique, and the Korea 
Science and Engineering Foundation. The Southeastern Universities Research Association 
(SURA) operates the Thomas Jefferson National Accelerator Facility for the United States 
Department of Energy under DOE contract DE-AC05-84ER40150 Modification No. M175. 



[1] W. Melnitchouk, M. Sargsian, and M. I. Strikman, Z. Phys. A359, 99 (1997), nucl-th/9609048. 

[2] L. L. Frankfurt and M. I. Strikman, Pliys. Kept. 76, 215 (1981). 

[3] F. E. Close, R. L. Jaffe, R. G. Roberts, and G. G. Ross, Phys. Rev. D31, 1004 (1985). 

[4] C. E. Carlson and K. E. Lassila, Phys. Rev. C51, 364 (1995), hep-ph/9401307. 

[5] C. E. Carlson, J. Hanlon, and K. E. Lassila, Phys. Rev. D63, 117301 (2001), hep-ph/9902281. 

[6] S. Simula, Phys. Lett. B387, 245 (1996), nucl-th/9605024. 

[7] J. J. Aubert et al. (European Muon), Phys. Lett. B123, 275 (1983). 

[8] M. M. Sargsian et al., J. Phys. G29, Rl (2003), nucl-th/0210025. 

[9] L. L. Frankfurt, M. M. Sargsian, and M. L Strikman, Phys. Rev. C56, 1124 (1997), nucl- 
th/9603018. 

[10] C. Ciofi degli Atti, L. P. Kaptari, and B. Z. Kopeliovich, Eur. Phys. J. A19, 145 (2004), 

nucl-th/0307052. 

[11] J. P. Berge et al., Phys. Rev. D18, 1367 (1978). 

[12] V. L Efremenko et al., Phys. Rev. D22, 2581 (1980). 

[13] A. Bodek et al., Phys. Rev. Lett. 50, 1431 (1983). 

[14] R. G. Arnold et al., Phys. Rev. Lett. 52, 727 (1984). 

[15] G. Bari et al. (BCDMS), Phys. Lett. B163, 282 (1985). 



32 



[16] D. M. Aide et al., Phys. Rev. Lett. 64, 2479 (1990). 

[17] S. D. Drell and T.-M. Yan, Phys. Rev. Lett. 25, 316 (1970). 

[18] S. Dieterich et al., Phys. Lett. B500, 47 (2001), nucl-ex/0011008. 

[19] S. Dieterich (Jefferson Lab E93-049 and MAMI Al), Nucl. Phys. A690, 231 (2001). 

[20] S. Strauch et al. (Jefferson Lab E93-049), Phys. Rev. Lett. 91, 052301 (2003), nucl-ex/0211022. 

[21] R. D. Ransome, Nucl. Phys. A699, 360 (2002). 

[22] P. G. Blunden and G. A. Miller, Phys. Rev. C54, 359 (1996), nucl-th/9602031. 

[23] D. H. Lu, A. W. Thomas, K. Tsushima, A. G. Wilhams, and K. Saito, Phys. Lett. B417, 217 

(1998), nucl-th/9706043. 

[24] T. Kafka et al. (E545 collaboration). Bull. Am. Phys. Soc. 28, 756 (1983). 

[25] E. Matsinos et al. (BEBC-WA59 collaboration), Z. Phys. C44, 79 (1989), and quotes therein. 

[26] B. A. Mecking et al. (CLAS), Nucl. lustrum. Meth. A503, 513 (2003). 

[27] A. V. Klimenko and S. E. Kuhn, CLAS Note 03-005, Jefferson Lab (2003). 

[28] M. Lacombe et al., Phys. Rev. C21, 861 (1980). 

[29] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C51, 38 (1995). 

[30] K. Abe et al. (E143), Phys. Rev. D58, 112003 (1998), hep-ph/9802357. 

[31] P. L. Anthony et al. (E155), Phys. Lett. B493, 19 (2000), hep-ph/0007248. 

[32] M. Arneodo et al. (NMC), Phys. Lett. B364, 107 (1995), hep-ph/9509406. 

[33] A. Bodek et al., Phys. Rev. D20, 1471 (1979). 

[34] L. H. Tao et al. (E140X), Z. Phys. C70, 387 (1996). 

[35] L. M. Stuart et al., Phys. Rev. D58, 032003 (1998), hep-ph/9612416. 

[36] K. Abe et al. (E143), Phys. Lett. B452, 194 (1999), hep-ex/9808028. 

[37] P. Bosted, Phys. Rev. C51, 409 (1995). 

[38] L. W. Mo and Y.-S. Tsai, Rev. Mod. Phys. 41, 205 (1969). 

[39] A. Afanasev, I. Akushevich, V. Burkert, and K. Joo, Phys. Rev. D 66, 074004 (2002), hep- 



ph/0208183, URL |http : //www . j lab . org/RC/exclurad/ 



[40] J.-M. Laget, Phys. Lett. B609, 49 (2005), and private communication. 

[41] C. Ciofi degh Atti and B. Z. Kopehovich, Eur. Phys. J. A17, 133 (2003). 

[42] C. Ciofi degli Atti, private communication. 

[43] C. Butuceanu (2005), to be published, Ph.D. thesis. College of William and Mary. 



33 



